Quantum Phase Transition in Hall Conductivity on an Anisotropic Kagome Lattice 
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We study the quantum Hall effect(QHE) on the Kagome lattice with anisotropy in one of the 
hopping integrals. We find a new type of QHE characterized by the quantization rules for Hall 
conductivity a xy = 2ne 2 /h and Landau Levels E(n) = ±df\/ (n + l/2)hBe (n is an integer), which 
is different from any known type. This phase evolves from the QHE phase with a xy = 4(n + l/2)e 2 jh 
and E(n) = ±VF\f2nhBe in the isotropic case, which is realized in a system with massless Dirac 
fermions (such as in graphene). The phase transition does not occur simultaneously in all Hall 
plateaus as usual but in sequence from low to high energies, with the increase of hopping anisotropy. 
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The quantum Hall effect (QHE) is a remarkable trans- 
port phenomena in condensed matter physics There 
are three kinds of integer QHE in the known materi- 
als. One is the conventional integer QHE occurring 
in two-dimensional (2D) semiconductor systems, where 
the successive filling of the Landau levels leads to an 
equidistant ladder of quantum Hall plateaus at inte- 
ger filling n — 0, ±1,±2, with a quantized value 
<j xy = 2ne 2 /h The second is the unconventional QHE 
observed in graphene, where charge carriers mimic the 
massless Dirac fermions, so that the Hall conductivity 
is half-integer quantized a xy = 4(n + l/2\e /h due to a 
Berry phase shift 7r at the Dirac points 0, y|, 0, [H, 0, 0, • 
The third occurs in bilayer graphene, where the charge 
carriers have a parabolic energy spectrum but are chiral 
with a Berry's phase 2ir. Therefore, the Hall conductiv- 
ity follows the same ladder as in conventional 2D electron 
gases, but the plateau at zero level is absent (il.[lol|. 

In this Letter, we demonstrate a new kind of integer 
QHE on the anisotropic Kagome lattice and the quantum 
phase transition relating it to the unconventional QHE 
for the massless Dirac fermions. The Kagome lattice has 
recently attracted considerable interest due to its higher 
degree of frustration. It is the line graph of the hon- 
eycomb structure in view of the graph theory [lit ]. The 
three-band electronic structure [Fig. 1] is composed of one 
flat band and two dispersive bands. The latter has the 
same form as that in graphene [3], and the two bands 
touch at two incquivalent Dirac points forming massless 
Dirac fermions. As a result, an unconventional QHE 
with a xy — 4(n + l/2)e 2 /h is realized on the isotropic 
Kagome lattice. Assuming one of the three hopping inte- 
grals, which is denoted by t%3, can take a different value 
from the two others, we find a quantum phase transition 
for the Hall conductivity from the unconventional form 
&xy = 4(n + l/2)e 2 /h to a xy = 2ne 2 /h. Though the 
latter phase shows the same Hall conductivity as that 
found in conventional 2D semiconductors, it has the fol- 
lowing non-trivial properties, i) The phase is character- 
ized by a new quantization rule of Landau levels (LL) 
with E(n) = VF\/{n + l/2)HBe, in contrast to E(n) — 



(n + \/2)Thjo c with uo c = eB/m for the free-fermion QHE 
systems, E(n) = VF\/2nhBe for the single-layer graphene 
and E(n) = y/ n(n — \)hw c for the bilayer graphene. ii) 
The quantum phase transition does not occur simultane- 
ously in all Hall plateaus as usual but in sequence from 
low to high energies with the increase of anisotropy. iii) 
The quantum phase transition occurs only in the case of 
^23 < *12 ( f i2 = hz). In the other case of i 23 > *12, 
the unconventional QHE realized in the isotropic system 
remains at least for the largest anisotropy we consid- 
ered here, namely £23 = 2ii 2 . This kind of quantum 
phase structures controlled by the anisotropy of the hop- 
ping parameters is also in stark contrast to that in the 
honeycomb lattice(graphene), where the unconventional 
QHE evolves into the conventional one for the strong £23 
(^23 > ^12) regime, while no phase transition occurs for 
the weak i 23 (*23 < £12) regime [12], |l3|, [lj, |l5j. There- 
fore, the quantum phase transition demonstrated here 
has no known analogues and presents an intriguing case 
for experimental studies on QHE. 

We start from the tight-binding model on a 2D metallic 
Kagome lattice, 



H 



(1) 



where Ci a {& irT ) annihilates(creates) an electron with spin 
0(0 =t,|) on site i and Uj is the hopping integral be- 
tween the nearest neighbors. Considering that there 
arc three sites in each unit cell [see Fig. Ufa)], we 
can write Eq.( Q]) in the momentum space as H = 

E kCT *L M ( k )*k ff . Where * ff (k) = (CklCT, Ck2cr, Ck3cr), 
and M(k) is a 3 x 3 matrix 



e 12 (k) Ei 3 (k)\ 
M(k) = ( e 2 i (k) e 23 (k) 
e 3 i (k) e 32 (k) / 



(2) 



with £12 (k) = £21 (k) = -2i i2 cos(k • Si), e 13 (k) = 
£31 (k) = -2i 3:L cos(k • S3) and £ 2 3(k) = £ 3 2(k) = 
— 2t 2 3COs(k • 82)- Si, 8 2 and S3 are the nearest-neighbor 
vectors, Si = (l/2)x, S 2 = (l/4)£ + (V3/4)y, S3 = 
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FIG. 1: (color online) (a) Lattice structure of the Kagome 
lattice. The region enclosed by the dash lines represents the 
Wigner-Seitz unit cell, ai and a.2 are the lattice unit vectors. 
The lattice sites in each unit cell are labeled by 1, 2, and 3. 
(b) The first Brillouin zone, bi and b2 are the reciprocal- 
lattice vectors, (c) The electronic dispersion for ti% = £23 
/ - = 1.0. The Dirac cones are located at K(Ki) and K'(K[) 
points as labeled in Fig.(b). 



— (l/4):i + (\/3/ 4)y. For the isotropic case (ti2 = £23 
£31 = t), the energy bands are, 



E±(k) = -t±Vl + /(k), E 3 {k) = 2t, 



(3) 



where /(k) = 8 cos(fc x /2) cos(fc x /4 + y/3k y /4) cos(fc x /4 - 
VSky/i), and the energy band structure is shown in 
Fig. [He). The two dispersive bands E + and E~ 
contact at the Dirac point K = (47r/3,0)[or K' — 
(27r/3, 2n/y/3)]. Near the Dirac point q = k— K, the elec- 
trons behave as Dirac fermions with the approximated 
dispersion, 



E ± =-t±v F \q\ + 0[(q/K)% 



(4) 



vp = 3v / 2t/4 is the Fermi velocity. Thus, except for an 
additional flat band £3, the two dispersive energy bands 
are similar to those in graphene. As a result, an un- 
conventional QHE as that found in graphene is expected 
on the isotropic Kagome lattice around the Fermi level 
li = —t (corresponding to the electron filling 1/3). 

In an uniform magnetic field applied perpendicular to 
the sample plane, the tight-binding Hamiltonian is, 



H 



E 



(tije 



' C i>r C ')Cr + 



H.C.). 



(5) 



The Hamiltonian Eq.(5) will be diagonalized numerically 
on a finite lattice with size N — 3 x L\ x L2 (factor three 
counts the three sites in each unit cell). The magnetic 
flux per triangular is chosen to be (j> = ^ A cty = with 

M an integer, then the total flux $ = ^jj- through the 




FIG. 2: (color online) (a)-(c) Hall conductivities for different 
t2t as a function of the Fermi energy and the corresponding 
density of states which is denoted by red vertical lines [(a) 
t 2 3 = 1.0. (b) t 23 = 0.4. (c) t 23 = 0.15]. The numbers 



on each plateau denote v in a 



xy = ve /h. The red arrows 
in (b) and (c) indicate the transition point from g 3 — 4 to 
g 3 = 2. (d) Phase diagram showing the gradual splitting of 
the Hall plateaus with the decrease of t23- The solid lines 
denote the original plateaus and the dashed lines the new 
plateaus induced by the anisotropy. 



lattice is taken to be an integer to satisfy the periodic 
boundary condition. Typically, N — 3 x 320 x 320 and 
<t> = tJzh are used in numerical calculations. After the 

f 2560 

diagonalization, the Hall conductivity a xy is calculated 
with the Kubo formula 



' xy 



A 



E 



\v x \P)Wy 



\v y \/3)(/3\v x \a) 



(6) 

where A = i2e 2 /SH with S the area of the system, v is the 
velocity, e a and Ep are the corresponding eigenvalues of 
the eigenstates |a) and \(3). In the following, the hopping 
integrals ii2 = £31 = 1.0 are used as the energy unit and 
t23 to acount for the anisotropy. 

In Fig. [21 the Hall conductivity <j xy — v^- and the 
electron density of states(DOS) are plotted for the dif- 
ferent £23- For the isotropic case with £23 = 1-0, the 
Hall plateaus satisfy the unconventional quantization 
rule v = (n + l/2)g s with a degeneracy factor g s = 4 
for each Landau levcl(LL) to count two spin components 
and two Dirac points. The 1/2 shift of the Hall plateaus 
is due to a nonzero Berry phase at the Dirac points [l6| , 
or can be simply explained as arising from the existence 
of a zero mode as seen from the DOS shown in Fig. 2 (a). 
This numerical result is consistent with the above analyt- 
ical calculation, showing that the QHE on the isotropic 
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Kagome lattice exhibits the same behavior as that in 
graphene. 

Decreasing i 2 3 to introduce the anisotropy in the hop- 
ping integral, one will find that the steps in Hall con- 
ductivity split at the mid-point gradually, as shown in 
Fig. 2(b) and Fig. 2(c). Concomitant with the splitting, 
a new Hall plateau emerges between every two Hall 
plateaus and the degeneracy factor g s changes from 4 
to 2. Basically, one may expect that, upon the intro- 
duction of the anisotropy, the rotational symmetry of 
the isotropic Kagome lattice will be broken, and con- 
sequently the degenerate energy levels will be separated. 
Indeed, at the energy levels where the step splitting in 
Hall conductivity occurs, the peak of DOS (denoted as 
red vertical lines in Fig. 2) is split into two adjacent peaks 
with half a previous height. However, two nontrivial 
characters exhibit here. One is that the splitting does not 
happen as a whole simultaneously, but gradually from 
low energies to high energies with the decrease of £23 (cor- 
responding to the enhancement of the anisotropy) . More- 
over, whenever a splitting occurs, a new Hall plateau 
emerges. But, the next splitting does not happen in 
succession with the further decrease of £23, instead the 
emerging plateau will grow firstly until it satisfies the 
new quantization rules for the Landau Level which will 
be addressed in the following. This process can also be 
clearly seen from the peak splitting in DOS. In this way, 
the QHE on the anisotropic Kagome lattice exhibits a 
sequent quantum phase transition from low to high ener- 
gies, as shown in the phase diagram presented in Fig. 2(d). 
Secondly, in the case of £23 > 1 (see Fig. [3fa)) ; the peaks 
in DOS does not split at any energy level in the reason- 
able parameter regime. As a result, at least for the largest 
anisotropy £23 = 2 considered here, the QHE shows the 
same behavior to the isotropic case. 

The quantum phase transition demonstrated above 
is in stark contrast to that on the honeycomb lat- 
tice (graphene) [HI, 13, Q, 15 1, where the unconven- 



tional QHE changes into the conventional QHE with 
(T xy = 2ne 2 /h in the strong £23 (^23 > 1) regime, while no 
phase transition occurs for the weak £23 (^23 < 1) regime. 
In addition, the phase transition on the anisotropic hon- 
eycomb lattice occurs symmetrically starting from both 
the low and high energies, and gradually approaches to 
the zero-energy level, so that it exhibits a particle-hole 
symmetry. 

Next, let us study the character manifested by the 
new phase emerging in the quantum phase transition. 
This is accessible easily in the limit £23 = 0, where the 
Kagome lattice is topologically equivalent to the lattice 
in Fig. [3fc) with only the nearest neighbor hoppings. 
The three energy bands on this lattice can be found 
as, E± — ±2t^cos 2 (k x /2) + cos 2 (fc y /2) and E 3 = 0, as 
shown in Fig. [3^d) . Around Ko = (n, 7r), the bands 
E± have the linear dispersion E±(q) = ±vp\q\ with 
Vp = \f2~t. Because the flat band £3 crosses the Kq point, 




FIG. 3: (color online) (a) and (b) are the Hall conductivity 
together with the DOS for £23 = 1 and £23 = 0, respectively, 
(c) The square lattice which is equivalent topologically to the 
Kagome lattice for £23 = 0. (d) The energy bands for the 
lattice in (c). 



the electrons do not behave as massless Dirac fermions. 
In this respect, the Hall conductivity shows a conven- 
tional behavior a xy = 2ne 2 /h. However, the energy of 
the LL follows a new quantization rule(Fig.3(b)). This 
can be obtained analytically by solving the low-energy 
Hamiltonian of the system, as given by 



Px Py 

H=-^ [Pi 

V2 \ pi 



(7) 



under a uniform magnetic field B by replacing the mo- 
mentum operator p by p + eA with A = B(—y, 0). The 
result turns out to be, 



E(n) = ±v F y/{n + l/2)HBe, n = 1,2,3, 



(8) 



This kind of distribution of the LL leads to a cor- 
responding distribution of the Hall plateaus, which 
is different from that in conventional semiconduc- 
tors [n + 1/2)/kj c , graphenes vpy2nhBe and bilayer 
graphenes yj n(n — 1) hoj c . 

To understand the property of the quantum phase 
transition demonstrated above, we show in Fig. 4 the evo- 
lution of the energy band with the hopping integral £23 
along the high symmetrical directions. First, we point 
out that the unconventional QHE is found numerically 
to be limited to a finite energy range (not shown here), 
which is the region from A to A' (points of the van Hove 
singularity) in the dispersions shown in Fig. 4. Outside 
this region, the QHE will exhibit a conventional behavior 
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FIG. 4: (color online) The energy bands for different £23- (a) 
The energy bands for weak £23 regime. Two typical energy 
bands are plotted along V — K — K' — T direction. The blue 
dash, and red dash-dot lines are corresponding to £23 = 0.4 
and i23 = 0.15 respectively, (b) The energy bands for strong 
£23 regime. The typical energy bands for £23 = 2(the red dash 
lines) are plotted along M — K\ — F — K2 — M' direction. In 
the two figures, the black solid lines are the energy bands for 
the isotropic case(f23 = 1). 



in Hall conductivity. For t 2 3 < 1 [Fig. 4(a)], the two Dirac 
points, which are at K and K' points in the isotropic case, 
approach each other along the K — K 1 direction with the 
decrease of £23- Interestingly, in this process, the energy 
band around the A point is suppressed and lifted upwards 
gradually, while that around the A 1 point is not changed. 
As a result, the low energy part will be excluded out of 
the unconventional regime. Therefore, the conventional 
Hall conductivity emerges at these energy levels. 

On the other hand, for £23 > 1 [Fig. 4(b)] the two Dirac 
points approach each other along the K\ — T — K[ di- 
rection [see Fig. 1(b) for illustration]. Different from the 
case of t23 < 1, the energy band around A now is shifted 
downwards with the increase of £23 > so that the energy 
region exhibiting the unconventional Hall conductivity 
is enlarged. Therefore, no quantum phase transition is 
observed in the energy region considered. 

Finally, let us give a few comments on the possible ex- 
perimental realization of the theoretical prediction elabo- 
rated here. The anisotropy of the hopping integral can be 
realized by the distortion of the lattice due to the mono- 
clinic distortion, such as in Cu(3)V207(OH) 2 • 2H2O 17j . 
or by the difference in orbital characters on the atomic 
sites in each unit cell due to the Jahn- Teller effect, such 
as in Cs2Cu3CeFi2 (l8| . On the other hand, the Kagome 
lattice has been proposed to be realized by implement- 
ing an optical lattice for ultra-cold atoms [H, [2(| ■ In 
this regard, the ability to conveniently control the phys- 
ical parameters in the system facilitates the realization 
of the anisotropy. It also worthwhile to point out that 
the 2D lattice in the limit £23 = as shown in Fig. 3(c) is 
different from the Kagome lattice, and may provide an 
interesting model system for experimental investigation 
on the special quantum dynamics demonstrated above. 

In summary, we have study the quantum Hall effect 



on the anisotropic Kagome lattice. The anisotropy is in- 
troduced by assuming one of the hopping integrals £23 
taking a different value. In the weak £23 (£23 < ^12) 
regime, we find a new type of QHE characterized by the 
quantization rules for Hall conductivity a xy — 2ne 2 /h 
and Landau Levels E(n) — ±UfV (n + l/2)hBe, which is 
different from the known types. This phase evolves from 
the unconventional QHE with a xy = 4(n + l/2)e 2 /h via 
a quantum phase transition, which occurs successively 
from low to high energies with the decrease of the hop- 
ping integral. This quantum phase transition is absent in 
the strong £23 regime. Possible experimental realization 
of the theoretical prediction is also discussed. 
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